10 PRINT : PRINT "EXAMPLE 3.11-14 EMPTYING OF RESERVOIR": PRINT
20 DEFINT I
30 READ V0, H, HF, G, L: DATA 8.02,35.52,20.,32.2,1000
40 PRINT : PRINT "V0,H,HF,G,L="; V0; H; HF; G; L: R = (HF - V0 ^ 2 / (2! * G)) / V0 ^ 2
50 C1 = G / L: C2 = R + 1! / (2! * G): C3 = -.01: DT = 5!: T = 0!: V = V0
60 OH = H: U11 = DT * C3 * V: U12 = DT * C1 * (H - C2 * V ^ 2)
70 U21 = DT * C3 * (V + U12): U22 = DT * C3 * (H + U11 - C2 * (V + U12) ^ 2)
80 H = H + .5 * (U11 + U21): V = V + .5 * (U12 + U22): T = T + DT: REM PRINT H; V; T
90 IF H < 20! THEN T = T - (20 - H) * DT / (OH - H) ELSE GOTO 60
100 PRINT "T="; : PRINT USING "####.##"; T; : PRINT " SECONDS"